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We study the interplay of natural supersymmetry and a supersymmetric solution to the dis¬ 
crepancy observed in measurements of the muon magnetic moment. The strongest constraints on 
the parameter space currently come from chargino searches in the trilepton channel and slepton 
searches in the dilepton channel at the LHC, and vast regions are currently allowed, especially at 
large tan jl. With light top squarks in the spectrum, as required from naturalness arguments, the 
situation changes dramatically; stop-assisted chargino and neutralino production via i hx^ and 
i txi are already probing the entire parameter space compatible with the muon magnetic moment 
at tan yd ~ 0(10), while upcoming stop searches will probe most of the parameter space at larger 
tan yd ~ 40. Direct Higgsino searches as well as compressed slepton and stop searches are crucial to 
close out corners of parameter space. We consider one such example: in the presence of light slep- 
tons and charginos as required to obtain appreciable contributions to the muon magnetic moment, 
compressed stops can dominantly undergo the following decay i biiy(i>£) facilitated by 

off-shell charginos. We find that the enhanced branching to leptons leads to a 5a mass reach (with 
3000 fb“^ of data at LHC 14) of ~ 350 GeV, with the mass difference between stops and the 
lightest neutralino being ~ 80 GeV. This will further close out a significant part of the parameter 
space compatible with naturalness and the muon magnetic moment. 


I. INTRODUCTION 


The anomalous magnetic moment of the muon, = 
{g — 2)^/2, is one of the most precisely measured quan¬ 
tities in particle physics [T]. It will soon be measured to 
even higher degrees of precision by the New g — 2 Ex¬ 
periment at Fermilab [2] and the J-PARC experiment in 
Japan [3]. Since its theoretical value within the Standard 
Model (SM) can be calculated to within sub-parts-per- 
million precision, a comparison between theory and ex¬ 
periment provides a powerful probe of new physics. In 
fact, there has been a long-standing discrepancy between 
the measured and SM values, denoted by Aa^, given by 

Aa^ = a^(exp) - a^(SM) = (26.1 ±8.0) x 10“^^ . (1) 

The above value uses [4] for contributions of the hadronic 
vacuum polarization, and [5] for the hadronic light-by¬ 
light contribution. The theoretical calculation by [6] 
gives the discrepancy to be Aa^ = (28.7 ± 8.0) x 10“^^. 
Either way, one has a ^ 3cr deviation from the SM value, 
providing a tantalizing possibility of new physics. 

Supersymmetry provides a particularly appealing 
framework to understand this anomaly [ 7 ]. The size of 
the discrepancy is comparable to the SM Electroweak 
contribution a^(EW) - (15.36 ± 0.1) x 10-^^ 
Parametrizing new physics contributions to Aa^ by a 
coupling constant anew and a new physics scale Anew in 
the loop, this implies that 


^new 


Aa^ ^ 
aEw, 


^new 

Anew = SUSYl ^ Aew • 


( 2 ) 


Much like the dimensional arguments leading to the so- 
called dark matter WIMP miracle, the fact that the muon 


g — 2 anomaly prefers new physics at the Electroweak 
scale makes supersymmetry a prime contender in its res¬ 
olution. The parameter space of the Minimal Supersym¬ 
metric Standard Model (MSSM) has been widely studied 
in this context m, m 

The supersymmetric particles that are most relevant 
for significant contributions to the ggj vertex (smuons, 
muon sneutrinos, charged and neutral Higgsinos and 
Winos, Binos) all have small direct production cross- 
section at the Large Hadron Collider (LHC). Hence, it 
is not a surprise that the first run has not yielded any 
evidence for a supersymmetric resolution of the muon 
g — 2 anomaly. The situation, however, is considerably 
more serious if one considers the main motivation for su¬ 
persymmetry - stabilizing the gauge hierarchy. One can 
then advance fine-tuning arguments to make the case 
that color-assisted production of these particles should 
already have probed a substantial part of the available 
parameter space. The quadratic divergence in the Higgs 
mass (coming dominantly from the top sector of the 
SM) is cancelled in supersymmetric extensions by loops 
of scalar partners of the top quark, the top squarks or 
stops. Stabilizing the Higgs mass with small fine-tuning 
thus requires light stops. Since the Electroweak symme¬ 
try breaking condition within the MSSM shows that the 
Higgs vacuum expectation value (vev) cannot be much 
larger than the Higgsino mass parameter g, natural Elec¬ 
troweak symmetry breaking also requires light Higgsinos. 


^ Although our focus will be supersymmetry, we point out that 
the muon magnetic moment has been studied in other contexts 
recently, for example in M- 
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This is the framework of natural supersymmetry m 

In this paper, we investigate the interplay between a 
natural supersymmetric spectrum and the possibility of 
addressing the muon g — 2 anomaly within supersym¬ 
metry, in the backdrop of constraints coming from the 
first run of the LHC. We should clarify that we do not 
consider a natural spectrum in the strict sense that all 
superpartners that don’t play a role in Electroweak sym¬ 
metry breaking are decoupled; that would exclude light 
sleptons, for example. Rather, we search for regions of 
parameter space that would be compatible with the ob¬ 
served discrepancy of the muon magnetic moment and 
the principal components of a natural spectrum - light 
Higgsinos and light stops. 

Contributions to the muon magnetic moment in the 
MSSM are enhanced in two cases: the chargino/sneutrino 
diagrams and the Bino/smuon diagrams. These contri¬ 
butions are constrained mainly by chargino searches at 
the LHC in trilepton final states (in the case that the 
chargino being directly produced is a charged Wino) and 
by slepton searches in dilepton final states. Broadly, we 
find that the combination of these searches is only able 
to rule out a modest part of the parameter space for 
tan^d ^ 0(10), while for larger tan^d ^ 40, the pa¬ 
rameter space compatible with Aa^ is essentially uncon¬ 
strained now, and will depend on the future performance 
in these two channels at 14 TeV. 

Taking into acount color-assisted production of these 
particles from light stops changes the situation dramati¬ 
cally. Stop searches in the i bxf channel are probing 
the entire parameter space in case of chargino/sneutrino 
contributions for tan^d ^ 0(10), and much of the pa¬ 
rameter space for tan /3 ^ 40. Similarly, stop searches in 
the t txi channel are probing much of the parameter 
space in case of Bino/smuon contributions, for a variety 
of choices of other parameters. The performance of the 
LHC in these two stop search channels at 14 TeV will 
probe much of the remaining parameter space compati¬ 
ble with natural supersymmetry and muon g — 2. This 
presents a sobering perspective on the prospects of ad¬ 
dressing both the gauge hierarchy and the discrepancy in 
muon g — 2 within minimal supersymmetry. 

If one is willing to allow the fine-tuned possibility that 
the mass difference between the stop and the lightest su¬ 
persymmetric particle (LSP) is less than the top mass 
(stop ’’three-body decay region”) or even less than the 
W mass (stop ’’four-body decay region”), then the con- 


^ On the other hand, it is the stop sector that provides the domi¬ 
nant radiative corrections to the Higgs quartic coupling; obtain¬ 
ing the measured value of 125 GeV for the Higgs mass within 
the MSSM thus requires stops above a TeV, unless the A-terms 
are tuned to maximal mixing. This tension between the fine- 
tuning associated with obtaining a 125 GeV Higgs, and at the 
same time obtaining natural Electroweak symmetry breaking, is 
the little hierarchy problem. Assuming that the little hierarchy 
problem can be solved through some extension of the MSSM, it 
is reasonable to expect a natural spectrum. 


straints on stops are considerably weaker and much of 
the parameter space available for the muon magnetic mo¬ 
ment opens up, as the bounds revert exclusively to the 
chargino and slepton searches. Our second set of main 
results is to probe this important caveat: the presence 
of light, compressed stops which could have evaded ob¬ 
servation. In the presence of light charginos and slep- 
tons, an important decay mode of the stop opens up: 
i biiy{i>i) biuxi- This decay can be 

competitive or even dominate over conventional 3-body 
or 4-body decays of stops through on-shell or off-shell 
W’s. We study this decay in the compressed regions, 
where 175 GeV > (m^ — m^o) > 85 GeV, and 85 GeV 

^ {mi — m^o) > 0 GeV. We find that the 14 TeV LHG 
can discover stops in this decay mode up to a mass of 
^ 350 GeV at a luminosity of 3000 fb“^, when the mass 
difference of the stop with the LSP is ^ 80 GeV. This will 
serve the function of further closing out the muon g — 2 
parameter space from the point of view of naturalness. 

The rest of the paper is structured as follows. In Sec¬ 
tion |n| we describe the features of the parameter space 
in the MSSM that give rise to appreciable contributions 
to the muon magnetic moment. In Section [ni| we sum¬ 
marize existing collider constraints on the superpartners 
that are relevant to our study. In Section EYl we give 
our first set of main results, showing the existing collider 
constraints and top squark probes in regions of parame¬ 
ter space that are compatible with the observed value of 
Aa^. In Section |V| we discuss the results from the pa¬ 
rameter space analysis. In Section |VTj we present our sec¬ 
ond set of main results, analysing the collider prospects 
of compressed stops decaying to the lightest neutralino 
via sleptons. We end with our Gonclusions. 


II. MUON ^ - 2 IN SUPERSYMMETRY 

The MSSM parameter space can easily accommodate 
the observed discrepancy of the muon magnetic moment. 
The one-loop contributions are well known and displayed 
in Fig. [2 The left panel shows the case where a smuon 
and a neutralino dominate the one-loop contribution. In 
this case, the charginos can be heavy. On the other hand, 
the right panel shows the case where the contribution is 
driven by a light chargino and a muon sneutrino. Since 
the LSP has to be a neutralino, this means that the spec¬ 
trum in this case has a light neutralino, in addition to a 
chargino and the muon sneutrino. 

Both classes of contributions can be probed by study¬ 
ing the space spanned by the following parameters: 

(3) 

where g is the Higgsino mass parameter, tan /3 is the ratio 
of the Higgs vev’s. Mi and M 2 are the Bino and Wino 
masses, and are smuon masses, and is the 
soft mass of the muon sneutrino. 
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The analytic expressions for these diagrams are given 
below mniHiii. 
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where ^ and g' are the gauge couplings of the SU (2) and 
U(l) SM groups, respectively, and the ] and 

are loop functions given by 


J 5 — 3(x + ^) + 2 InX 

\ (x —1)^(^ —1)^ x — y {x — 1)^ 
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= xy 
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In the above, the reduced forms of [18] have been 
used. Also, A~±~ denotes contributions from 

Xl 

chargino/sneutrino diagrams, while A^^~, A^^~, A^^~, and 

A^^~ denote the contributions from neutralino/smuon di¬ 
agrams. 

Several features of the above equations should be 
noted: 


• Larger tan (3 values will result in the most conser¬ 
vative collider bounds on the superpartners. This 
is because large tan/3 enhances the supersymmet¬ 
ric contribution to the magnetic moment. This 
enables the superpartners participating in the di¬ 
agrams to have higher masses, where the collider 
bounds are weaker. In our study, we will generally 
use tan = 40. 


light, while the other can be heavy. This has some 
implications for collider bounds, which are typically 
more severe for left-handed sleptons. 

• In the limit of heavy charginos, the neutralino- 
smuon contribution of Eq. dominates. Again, 
either the left or right-handed smuon has to be 
light; depending on how heavy the charged Wino 
and Higgsino are, they may both have to be light. 


In the next Section, we go on to discuss the relevant 
collider bounds coming from the LHC, before applying 
the bounds on the MSSM parameter space compatible 


with Aa^ in Section IV 



"I FIG. 1: Feynman diagrams contributing to muon g — 2 
I at one-loop in the MSSM. The left diagram shows the 
case involving a neutralino and a smuon, while the right 
diagram shows the case involving a chargino and a 
muon sneutrino. 


III. COLLIDER CONSTRAINTS 

In this section, we discuss the limits on supersymmet¬ 
ric spectra obtained from LEP and LHC searches. This 
will enable us to conveniently refer back to the relevant 
search when we discuss the parameter space in the next 
section. We will concentrate on the particles relevant for 
us, namely Xi-, f and Xi- 


LEP bounds: 


• The chargino-sneutrino contributions are sup¬ 
pressed for large /i, i.e. heavy Higgsinos. In fact, we 
will find later that the optimal contributions from 
these diagrams to Aa^ come from Higgsinos with 
mass ^ 0(500) GeV, which is precisely the region 
preferred by small Electroweak fine-tuning argu¬ 
ments. For heavier Higgsinos, either the charged 
Wino has to be light, or the neutralino-smuon con¬ 
tributions of Eq. I^have to dominate. 

• For the neutralino/smuon contributions, either the 
left-handed or the right-handed smuon has to be 


LEP results [19] constrain slepton and chargino 
masses: 

• sleptons: > 100.0 GeV where i = e, y; , 

• lightest chargino : m~± > 103.5 GeV. 

Xl 

LHC bounds: 

The search for gluinos {g) and first two generation 
squarks (g) at the Large Hadron Gollider (LHG) has so 
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far yielded null results. The exclusion limits on squark 
and gluino masses, when they are comparable, are ap¬ 
proximately 1.5 TeV at 95% CL with 20 of inte¬ 
grated luminosity [2QH23] . We will assume, in the re¬ 
mainder of the paper, that these states are decoupled 
from the low energy physics. Therefore, in studying the 
constraints on the superpartners relevant for our study, 
we will not consider the limits arising from gluino/squark 
production followed by cascade decay to Xi^ ^5 oi* xf • As 
discussed in the Introduction, cascade decays of the top 
squark will be important for us. 

(I) Chargino and neutralino masses - Limits on the 
lightest charginos coming from direct chargino produc¬ 
tion or production from top squarks are as follows. 

(a) We first consider the case where the Bino is the 

LSP, and there is no light slepton in the spec¬ 
trum, i.e., > m~±,m^o, where the neutral 

and charged Winos constitute X 2 xf 5 respec¬ 
tively. For direct charged Wino production, the 
signatures are 01 final states, 

with model dependent branching fractions. Tak¬ 
ing into account both ATLAS and CMS results, 
one can rule out charged Wino masses below 350 
GeV for Xi ^ 140 GeV [21]. 

(b) We next consider the case where the Bino is the 

LSP, and there is a light slepton in the spec¬ 
trum, i.e., < m~±,m^o, where the neutral 

and charged Winos constitute X 2 nnd yf, respec¬ 
tively. For direct charged Wino production with 
a lighter slepton, the exclusion bounds are much 
stronger, coming from trilepton searches mus]. 
For Xi ^ 300 GeV, xt is excluded below 720 
GeV [m [ 25 ] assuming sleptons intermediate to 
xt and xt Compressed regions of paramemter 
space with m~± —m^o ^ 30 GeV are still allowed 

Xi 

due to insufficient pr- We will use Figures 7(a), 
7(6), 7(c), and 7{d) of [24] for our work. 

(c) For cases where the neutral and charged Winos 
themselves comprise the LSP and next-to-lightest 
supersymmetric particles, the bounds are much 
weaker and it is useful to think in terms of high 
luminosity studies at 14 TeV. Taking a back¬ 
ground systematic uncertainty of 1%(5%), with 
3000 fb“^ of data, the 14 TeV LHG is sensitive 
to Winos of 240 GeV (125 GeV) and Higgsinos 
of 125 GeV (55 GeV) [26U28] in a pure cut and 
count analysis. Shape analysis can improve the 
mass reach significantly m- 

(d) We now consider bounds coming from direct stop 
production, in spectra with mi > m ~±, m^o . For 
charginos produced on-shell from top squarks, 
there are very stringent limits depending on 
the relative mass separations of xf 5 ^^4 i. 
For large mass separations between the lightest 


chargino and the LSP, xt — ‘^xt direct produc¬ 
tion of stops followed by the decay i bxt rules 
out stops up to ^ 500 GeV with xt ^ 400 GeV 

[29] , [30]. 

(e) For mi > rn ~±, m^o and for compressed 
chargino-neutralino system with rn ~± — m^o ^ 5 
GeV, direct production of stops followed by the 
decay i bxt rules out stop masses up to mi ^ 
600 GeV for m^o ^ 200 GeV, unless chargino- 
stop is also compressed, with a mass separation 
of ^ 40 GeV [ 29 ]. For compressed chargino- 
neutralino system but with m^o ^ 200 — 250 
GeV, top squarks are ruled out up to 600 GeV 
except for a window 250 — 350 GeV. 

(II) Slepton masses - The limits on sleptons are particu¬ 
larly strong when the initial particle in the decay chain 
is a squark/gluino or a stop or chargino, as discussed 
previously. However, for direct production of sleptons, 
the limits are much weaker due to small production 
cross section. 

(a) The direct production of sleptons has been 

probed both by dilepton searches at ATLAS m 
and GMS [36]. The decay chain is pp It 
/+/“x5x?, with Br{l fmol 

states containing opposite-sign same-flavor non¬ 
resonant dileptons and missing transverse energy 
(pr)- The mass separation m^ — m^o is an im¬ 
portant factor in the exclusion plots. The mass 
reach with m^o = 0 GeV is m^ ^ 280 GeV and 
330 GeV for GMS and ATLAS respectively. We 
will use Figures 8(a), 8(6), and 8(c) of [31] for 
our analysis. 

(b) Gompressed spectra with smaller mass separation 
between i and x5 have eluded these probes. 
For mass difference of £ and Xi of 3-15 GeV, a 
recent study with VBF tagged jets and pY has 
found that the 14 TeV LHG may be sensitive to 
sleptons up to mass m^=115-135 GeV, with 3000 
fb“^ of data [32]. Gompressed slepton searches 
have been studied by several other groups as well 

m- 

(HI) Stop masses - The bounds on the mass of the lightest 
top squark (t) are less stringent than those on gluino 
or first two generation squark masses, due to smaller 
production cross section. Exclusion limits in the m^- 
m^o plane have been obtained for a variety of decay 
modes: 

(a) The simplest scenario for i studies is to con¬ 
sider the direct QGD production of i pairs with 
100% branching i txt Exclusion limits in the 
mi-m^o plane have been obtained in this decay 
mode (monolepton search [29], dilepton search 

[30] ). The limits depend on the mass differ¬ 
ence mi — (m^o + mt). Eor m^o < 140 GeV, 
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stop masses up to ^ 650 GeV have been 
ruled out, except for highly compressed scenarios 
^ 20 GeV. For m^o < 140 — 200 
GeV, the same limits of stop masses apply, except 
mi in the interval 320 - 420 GeV also opens up. 

(b) Three body stop decays into hW'x^ have also been 

bounded between the limits + mw + 

m^o and < mt ^ m^o. These constraints are 
weaker due to smaller pY and softer 5-jets. The 
limits extend up to ~ 300 GeV and m^o ^ 
120 GeV. While most of the parameter space in 
this region has been ruled out, regions close to the 
boundary -1 -mw+ are still allowed. 

(c) Four body decays of stops satisfying + 

mw + have been studied, and the current 
constraints reach m^ ^ 240 GeV and m^o ^160 


A. Reinterpreting Slepton and Chargino Bounds 
for Higgsino LSP 


In the LHC constraints for charginos (trilepton chan¬ 
nel) and sleptons (dilepton channel) described above, the 
LSP has been assumed to be a pure Bino. In describ¬ 
ing scenarios compatible with natural supersymmetry, we 
will often be interested in scenarios where the LSP is a 
Higgsino [34]. The bounds should therefore be reinter¬ 
preted for that case. 

The main difference (for left-handed sleptons) will 
come from the fact that now the sneutrino can have non¬ 
zero branching to charged leptons, while the slepton can 
have non-zero branching to neutrinos. The branchings 
depend on the Bino or Wino components of Xi^X 2 ^ 
xf, and are given by [35] 


r 4 ^ ixU) = Cl (rnz 
r(^' ^ 

viXi,2) = Cl {mz 

where 




- rnz- 
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+ (10) 
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IGtt 4s^c^ 


mn 


{m?p — m^y 


mn 


(13) 


(14) 


In the above, mj> and mj) being the the parent slep¬ 
ton mass and daughter neutralino/chargino mass, re¬ 
spectively. In the Bino LSP case for which the LHC 


bounds are given, only the slepton production is consid¬ 
ered, while the sneutrino branches completely into miss¬ 
ing energy. Now, both slepton and sneutrino produc¬ 
tions have to be considered, as well as their respective 
branchings to charged leptons. We will mainly consider 
Ml ^ M 2 , and ^ ^ 0.1—0.5. In that case, the combined 
production cross section times branching to charged lep¬ 
tons CF X Br reduces approximately to the Bino LSP case. 

We will thus use the same exclusion bounds, for Con¬ 
straints (lb) and {Ha), as those given by the collabora¬ 
tions. We mention, however, that for special choices of 
Ml, M 2 , and their relative signs, the slepton mass reach 
may be reduced to ^ 230 GeV or enhanced to ^ 470 GeV 

Ea. 

The story for the trilepton searches for charginos is 
similar. The main point is that the trilepton searches 
will target the production of the second lightest chargino 
xf, which is the charged Wino. This xf will then decay 
to a slepton or a sneutrino, which will then decay down to 
three possible states (the neutral Higgsino, and the two 
charged Higgsinos), all of which behave effectively as the 
LSP since the mass difference between the charged and 
neutral Higgsinos is small. The relevant decay chains are 
xt ,’^xf, where in the 

last step we have shown all the possible final states. 

What is relevant for collider bounds is that the ef¬ 
fective branching of the charged Wino to the LSP will 
change relative to the canonical case where the LSP is a 
Bino and the Higgsinos are heavy, which is the case for 
which the Gollaborations have given their results. We 
have computed the new branchings in the region of pa¬ 
rameter space relevant for us following Eq. 9-12, and 
found that the net branching is similar to the canonical 
case. We have thus used the bounds from the Gollabo¬ 
rations, which are actually a little conservative for our 
case. 


IV. Aa^ AND A NATURAL SPECTRUM 

In this section, we describe the parameter space regions 
which can account for the observed value of Aa^, in the 
presence and absence of a top squark in the low energy 
spectrum. 

In Fig. we show muon g — 2 constraints for a spec¬ 
trum containing light Higgsinos, sleptons, and stops. On 
the left panel, we show the case of tan^d = 10, with 
M 2 = 2(1. is kept heavy and decoupled. The verti¬ 
cal axis plots the Higgsino mass parameter g, while the 
horizontal axis plots the mass of the left-handed sleptons 
m ^^. We note that /i sets the scale for both the lightest 
neutralino LSP, as well as the lightest chargino, which 
have a mass separation of 0(5) GeV. The hatched region 
is ruled out from LEP bounds. The region between the 
red solid lines shows the region of parameter space that 
is within 2cr of the observed value of Aa^ . The blue solid 
contour shows the limits on the {m^^, m^o ) plane com¬ 
ing from the ATLAS dilepton study £l which 
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FIG. 2: The interplay of the muon g — 2 anomaly and a natural supersymmetric spectrum, plotted on the (/r, 

plane. M 2 = 2/i, while Ir is kept heavy and decoupled. On the left panel, tan^d = 10, while on the right panel, 
tan/3 = 40. The area between the red solid curves is consistent within 2cr of the observed value of Aa^. The blue 
contour shows the ATLAS dilepton bounds coming from Ir ^ which is Constraint (//a). The green region 
shows the part of Aa^ - compatible space that is constrained by the ATLAS trilepton searches coming from 
xt 2 55 5 which is Constraint (lb) adapted to the case of Higgsino LSP. The 

hatched region is ruled out by LEP constraints on charginos and sleptons. To the right of the black solid line 
marked (3), the LSP is a neutralino. The parameter space below the dotted lines is being probed by current stop 
limits. The black dotted line marked (1) marks the chargino/neutralino mass constraint for the case of stop-assisted 
production, which is Constraint (Je). The black dotted line marked (2) shows the future mass reach for direct 
chargino/neutralino production in the absence of light stops, with 3000 fb“^ of data at LHC 14, which is Constraint 

{lie). All masses are in GeV. 


is Constraint {I I a). The green region shows the part of 
Aamu - compatible space that is constrained by the AT¬ 
LAS trilepton searches coming from charged Wino decay, 
which is Constraint (lb). Following our discussion in Sec. 
iHAl we use the trilepton bounds on charged Winos from 
ATLAS, even in the case where the Higgsinos and not the 
Binos constitute the LSP. We note that this is a conserva¬ 
tive assumption for the region in parameter space shown 
in Fig. 

The left panel of Fig. has two dotted lines 

below which the parameter space is being probed 
by current stop searches: (1) This line shows the 
chargino/neutralino mass constraint in the case of 
stop-assisted production, which is Constraint (Je); (2) 
This line shows the future mass reach for direct 
chargino/neutralino production in the absence of light 
stops, which is Constraint (Jc). In addition, the solid 
black line marked (3) demarcates the regon where the 
lightest neutralino is the LSP, as required from dark mat¬ 
ter considerations. 

It is clear that much of the parameter space required 
for obtaining the correct value of Aa^ is already ruled 
out by the ATLAS dilepton and trilepton studies, while 
the entire parameter space is in fact being probed in the 
case of stop-assisted production (except the case of com¬ 
pressed stops). 


The right panel of Fig. shows the same parameter 
space, with tan/3 = 40 and M 2 = 2/i. £r is decoupled 
as before. In contrast to the left panel, we now see that 
the enhancement of Aa^ due to larger tan^d results in 
an allowed region between the red solid lines where the 
LSP and slepton masses are higher. The ATLAS dilepton 
constraints in blue are now irrelevant, as are the trilepton 
charged Wino constraints that have not been displayed. 
There is a small allowed region between m^o = 280 GeV 
and m^o = 400 GeV which lies above the current bounds 
coming from a light stop. This region will be explored at 
Run II of the LHC. 

Fig|^shows the same parameter space as Fig. except 
now M 2 is also decoupled from the low-energy spectrum. 
The left panel shows the case of tan/3 = 40, M 2 = 1 
TeV, and = 3 TeV being heavy and decoupled. The 
right panel shows the case of tan /3 = 40, M 2 = 1 TeV, 
and . The regions between the solid red lines 

are the allowed space for Aa^. In the left panel, the AT¬ 
LAS dilepton bounds rule out a significant portion of the 
parameter space. The trilepton chargino searches are ir¬ 
relevant since the charged Wino is too heavy. However, 
current bounds on stop-assisted production are probing 
the entire space. The presence of a light right-handed 
slepton in the right panel allows for larger values of m^o 
and . The dilepton constraints are mostly irrelevant 
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FIG. 3: The interplay of the muon g — 2 anomaly and a natural supersymmetric spectrum, plotted on the 

plane. The fixed parameters are M 2 = 1 TeV, tan (3 = 40. On the left panel, = 3 TeV is heavy and 
decoupled. On the right panel, The legend is the same as Fig. The hatched region is ruled out by 

LEP constraints on charginos and sleptons. All masses are in GeV. 


here, although current bounds on stop-assisted produc¬ 
tion still probes most of the parameter space. The re¬ 
maining space above = 280 GeV will be probed in 
Run II of the LHG. 

We next give results for the case where fi is large, but 
there is a light stop in the spectrum. A representative 
example of the constraints is shown in the left panel 
of Fig. We have taken tan^d = 40, /i = 2 M 2 , and 
Ml = ^M 2 . Lower tan^d will make the collider bounds 
more stringent, since the particles will be pushed to lower 
masses to account for Aa^. The LEP excluded regions 
are hatched. The corridor between the solid red lines 
indicates the region that is compatible within 2cr of the 
measured value of Aa^. The green dot-dashed line in¬ 
dicates the full region ruled out from trilepton searches 
with xf which is Constraint 

(/6). The blue solid line shows the limits from the AT¬ 
LAS dilepton study, which is Constraint (//a) as in the 
previous figures. Eor > m~±, the limits are weaker 
and mo del-dependent. To the right of the solid black 
line marked (4) is the region where the LSP is a neu- 
tralino, while to the right of the solid black line marked 
(5) is the region where • Erom top to bottom, 

the dotted black lines give chargino mass bounds for the 
various scenarios described in the text: (1) top squark 
m 2 > m~± + m^, with decay chain i 6%^, ruling out 
m~± < 400 GeV, as given by Constraint (Id); (2) de¬ 
coupled top squark and heavy slepton, with xf Whxi 
or xf WZxi^ ruling out m~± < 350 GeV for m^o ^ 
140 GeV, which is Constraint (/a); (3) top squarks and 
Ml decoupled, with the neutral Wino being the LSP, 
which is Constraint (/c). 

On the right panel of Eig. we show the constraints on 
the (Ml, plane for the case of M 2 = 2Mi and /i = 2 


TeV, with tan/3 = 40. The right-handed sleptons are 
kept fixed at = 1.5 x The region between the 

red solid lines is allowed by Aa^. The blue solid contour 
shows the ATLAS dilepton bound as before, while the 
green region is ruled out by ATLAS chargino searches. 
The black dotted line marked (1) shows the limits on 
m^o coming from stop decay i ty?, which is Constraint 

ruling out the parameter space that is below it. 
To the right of the solid line marked (2) is the region 
where the LSP is a neutralino. 

In Eig. we continue to show the parameter space in 
cases where /i, M 2 are heavy. On the left panel, we plot 
the case of tan /3 = 10, with g = M 2 = 5 TeV. The blue 
solid contour shows the ATLAS dilepton bound, which 
is Constraint {I I a). The black dotted line marked (1) 
shows the limits on m^o coming from stop decay i 
txi which is Constraint {Ilia), ruling out the parameter 
space that is below it. To the right of the solid line 
marked (2) is the region where the LSP is a neutralino. 
The right panel shows the same plane, with tan/3 = 40, 
/i = 3.3 TeV, and M 2 = 750 GeV. The value of M 2 was 
chosen to correspond to just above the limit where the 
ATLAS trilepton search, which is Constraint (lb), begins 
to be important. 


V. PROSPECTS: CURRENT AND FUTURE 

In the previous sections, we have studied the parameter 
space for the muon g — 2 anomaly in the MSSM. The 
following conclusions can be drawn: 

• If the low energy spectrum consists only of Hig- 
gsinos and sleptons, then the current LHC slepton 
searches rule out a part of the parameter space for 
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FIG. 4: The interplay of the muon g — 2 anomaly and a natural supersymmetric spectrum, plotted on the 

plane (left panel) and plane (right panel), with g kept large. Left panel: tan^d = 40, g = 2 M 2 , and 

Ml = I M 2 . The dot-dashed green contour is the entire region ruled out by the ATLAS trilepton searches with 
xf which is Constraint (/6), while the solid blue contour is ruled out by ATLAS dilepton 

searches, which is Constraint (//a). To the right of the solid black line marked (4), the LSP is a Bino. To the right 
of the solid black line marked (5), the slepton is heavier than the charged Wino. The area between the red curves 
gives the region compatible within 2cr of the observed value of Aa^. The regions below the black dotted lines are 
being probed by various searches (see text). Right panel: M 2 = 2Mi, g = 2 TeV, and tan/3 = 40. The 
right-handed sleptons are kept fixed at = 1.5 x The blue and red contours denote the same regions as the 
left panel. To the right of the solid black line marked (2), the LSP is a Bino. The dotted black line marked (1) 
denotes the limits on m^o coming from i txi^ which is Constraint {Ilia). The hatched region is ruled out by 
LEP constraints on charginos and sleptons. All masses are in GeV. 






FIG. 5: The interplay of the muon g — 2 anomaly and a natural supersymmetric spectrum, plotted on the (Mi,m^^) 
plane (left panel) and {Mi, mi ) plane (right panel), with g kept large. On the left panel, tan^d = 10, with 
g = M 2 = 5 TeV and = rni^. On the right panel, tan^d = 40, /i = 3.3 TeV, M 2 = 750 GeV, and Mi ^ 

The red and blue contours are the same as in previous figures. To the right of the line marked (2) is the region 
where the LSP is a Bino. Below the dotted black line marked (1), the space is being probed by stop-assisted Xi 
production, with limits on m^o coming from stop decay i ty?, which is Constraint {Ilia). The hatched region is 

ruled out by LEP constraints on sleptons. All masses are in GeV. 


tan^d ^ 0(10), while for larger tan/3 ^ 0(40), the 


parameter space is completely unconstrained. 
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• If the low energy spectrum consists of light stops in 
addition to Higgsinos and sleptons, then large por¬ 
tions of parameter space consistent with Aa^ can 
be probed by stop searches. The entire parameter 
space consistent with Aa^ for tan/3 = 0(10) has 
been probed; significant parts have been probed for 
larger tan ^ 40, as we saw for various choices of 
M 2 and . The high luminosity LHC will essen¬ 
tially probe this region. 

• In the case of large fi and light Wino (Fig. but 
with no light stops, large parts of parameter space 
corresponding to regions where > m ~± are al¬ 
lowed, primarily due to the weakness of LHC con¬ 
straints in that region. Taking into account stop 
decays, however, most of the parameter space has 
been probed. 


(Je). This constraint is especially relevant when 
the chargino/sneutrino contribution to muon g — 2 
is dominant. This decay mode has been studied by 
several authors [37]. 

2. i ^ txi, with mi ^ which is impor¬ 

tant for Fig. right panel, and Fig. where 
the Bino/smuon contribution to the magnetic mo¬ 
ment is dominant. This compressed region has been 
studied by several authors [38|, [39] . 

3. The conventional 3-body and 4-body stop decay 
regions: 175 GeV > (m^ — m^o) ^ 85 GeV, and 

85 GeV ^ (m^ —m^o) ^ 0 GeV, respectively. In 
these regions, the 6 -jet pt is typically too small to 
allow efficient 6 -tagging, in addition to small 
The 3-body region has been explored by [39] . 


• For the case of dominant Bino/smuon contributions 
to Aa^ (Fig. 1^, stops probe a significant part of 
the parameter space. 

• We note that for the case where the Bino/smuon 
contribution is dominant and both charginos are 
heavy, there is no constraint at all from the LHG 
currently in the large tan /3 ^ 40 regime, as is clear 
from the right panel of Fig. 5. The only constraint 
comes from a light stop, since the area below the 
line marked ( 1 ) is being probed currently. 

Gurrent limits on the Electroweak sector coming from 
the GMS and ATLAS collaborations can only constrain a 
small part of the avaiable space for Aa^. Taking into ac¬ 
count top squark decays, however, the situation changes 
dramatically. It is clear that a natural supersymmetric 
spectrum can tell us a lot about the available parameter 
space for muon ^ — 2 , and future probes of the top squark 
are important in this regard. 


A. Compressed Stops 

The limits on charginos and neutralinos coming from 
stop decay break down if the stop is compressed, i.e., if 
the total mass of the decay products of the stop is approx¬ 
imately equal to the stop mass. These are the scenarios 
where there is little in the final state. While top 
squark searches at the next run of the LHG will sweep 
out larger and larger parts of the parameter space for 
Aa^, it is important to note that compressed stops pro¬ 
vide scenarios that could have evaded detection, and are 
compatible with natural supersymmetry and muon g — 2. 

The remainder of the paper is devoted to studying such 
a scenario. The possible compression scenarios for stops 
are: 


1. i ^ bxf^ with mi ^ which is important 

for the natural supersymmetric spectrum in Fig. [^ 
and Fig. and described in Constraint (Id) and 


While all of these scenarios are important, we will es¬ 
pecially be interested in scenarios that are affected by 
the presence of light charginos and sleptons, which are 
required for a large contribution to the muon magnetic 
moment. In particular, the decay of stops to off-shell 
charginos 

i ^ biuxi (15) 

can be competitive with decays in the conventional 3- 
body and 4-body regions, which occur through on-shell 
and off-shell W’s. In the remainder of the paper, we will 
concentrate on probing the decay in Eq. 


VI. STOP TO SLEPTON DECAY 


Our benchmark scenario will be chosen such that the 


decay in Eq. 15 is competitive. The benchmark points 
are chosen with the criterion that it is compatible with 
all existing searches and gives the required value of (^f — 
2)^, following the right panel of Eig. The benchmark 
spectrum is given in Table The case of AM = (m^ — 
m^o) = 80 GeV, which will turn out to be our best-case 
scanerio, is shown. 


TABLE I: Relevant SUSY masses (in GeV) at 
benchmark point, tan^d = 36, and Aa^ = 25.9 x 10“^^. 
We keep X2 ^ xf = 748 GeV, h = 730 GeV, 61 = 760 
GeV, g = 3.3 TeV. All other colored and non-colored 
states are heavy. 


Particle Mass (GeV) Br 


ii 

348 

38% {Mvt), 62% {Wvt) 

i 

302 

100% {ix\) 

i> 

290 

100% (i/x?) 


594 

100% {lx\) 

Xl 

270 
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The challenges of probing the decay mode in Eq. 
with the spectrum of Table [T| are symptomatic of com¬ 
pressed stop searches in general: the lack of ^ (near 
the compressed limit, the Xi provides little transverse 
missing energy) and the softness of the final state 6-jets, 
which makes 6-tagging difficult. Indeed, the current con¬ 
straints on stops are feeble in the compressed regions, 
where 175 GeV > {m^ — m^o) ^ 85 GeV with dominant 

decay mode t bWxi^ GeV ^ {m^ — m^o) ^ 

0 GeV, with dominant decay mode i bW^ *)xi- The 
decay chain in Eq. [^has the advantage of full branching 
to leptons. 


A. Results and Analysis 

Although the current run of the LHG is at 13 TeV, we 
carry out our analysis at 14 TeV, since our sensitivity is 
only at very high luminosity. Inclusive tt* + jets samples 
are generated with i masses in the range of 200—400 GeV, 
choosing several different values of Xi ^rid hence AM. 
The left handed sleptons are kept at a mass diference 
^ 30 GeV above the Xv The Xi is mostly Bino. 

The spectrum is generated with SuSpect m . Signal 
and background samples are generated with MADGRAPH5 
im followed by the parton showering and hadroniza- 
tion with PYTHIA [42] and the detector simulation using 
DELPHES [43]. The study is performed in the 2l-\-2b-\-pT 
final state. The simulated backgrounds are tt + (0 — 2)jf, 
V -|- (1 — 3)j, and t + (1 — 3)j, where V stands for W 
and Z. We have used backgrounds from the Snowmass 
Energy Erontier Simulations [44] . 

The following selections are applied: 

(1) Preselection: Two isolated leptons with pt > 10 
GeV and two 6-jets with pt > 30 GeV in \r]\ < 2.5 are 
required. The 6-tagging efficiency is taken to be the de¬ 
fault DELPHES value. 

(2) Lepton pt’ The leading lepton is required to have 
Pt i^i) < 50 GeV, while the next-to-leading lepton is 
required to have prih) < 25 GeV. These selections are 
based on Eig. [^ left panel. 

(3) Events are required to satisfy > 0.25. 

Here, Mg// denotes the scalar-summed transverse mo¬ 
menta of jets and pr- This selection is based on Eig. 
right panel. 

(4) 6-jet Pt’ The leading 6-jet is required to have 
Pt(6i) < 60 GeV, while the next-to-leading 6-jet is re¬ 
quired to have Pt(^ 2) < 40 GeV. These selections are 
based on Eig. [^ left panel. 

(5) Events are required to satisfy mT 2 < 290 GeV, 
for an LSP selection of 270 GeV. This selection is based 
on the right panel of Eig. We have used the following 
definition of mT 2 

^T2 (PT^i , PT^2 = 

mm (max (ruT (pt^i, Ptp^), ^t(pt^ 2 ^ Pt, 2 ^^))) (1^) 


with the minimization carried out over all decomposi¬ 
tions of pY subject to + Pt, 2^^ = Pr- Here, ttit 

denotes the transverse mass, while ’Prii and denote 
the transverse momentum vectors of the selected leptons. 
pY is the missing momentum vector. In the calculation, 
the mass of missing momentum is chosen to be the mass 
of each signal point. 

In Eig. [^ left panel, we show the pt distribution 
(normalized to unity) of the leading and next-to-leading 
leptons for signal (red unshaded histogram) and back¬ 
ground (grey shaded histogram) after preselection cuts. 
The right panel shows the distribution (normalized to 
unity) of pY/Meff for signal (red unshaded histogram) 
and background (grey shaded histogram) after lepton se¬ 
lections. 

In Eig.[^ left panel, we show the Pt distribution (nor¬ 
malized to unity) of the leading 6-jet for signal (red 
unshaded histogram) and background (grey shaded his¬ 
togram) after the ^/Me// selection. Eor the signal, the 
leading lepton and 6-jet transverse momenta are approx¬ 
imately determined by the mass difference between {£, z>) 
and Xi^ and i and (^, z>), respectively, and are softer than 
the background values. The right panel shows the distri¬ 
bution (normalized to unity) mT 2 (right panel) for signal 
(red unshaded histogram) and background (grey shaded 
histogram), after all other selections. 

The cut flow table with corresponding cross-sections 
at each stage are shown in Table [H] for the benchmark 
point. We only display tt + j background, although we 
have taken into account all the backgrounds mentioned 
previously, which are subdominant. After all the cuts, 
the cross section of t -h (1 — 3) jets is 0.023 fb; the cross 
section of V -h (1 — 3) jets is 8.4 fb. 

TABLE H: Summary of the effective cross-sections (fb) 
for different benchmark signal points as well as the total 

background at LHG 14. Masses and momenta are in 
GeV. 


Selection 

Signal 

tt -h j 

Preselection 

7.4 

16888 

Lepton Pt 

4.2 

3251 

pr/Meff > 0.25 

3.1 

1697 

6-jet PT 

0.9 

87.5 

mT2 < 290 

0.74 

66.9 


The significance S/S + H, where S and B are the sig¬ 
nal and background rates, respectively, is plotted in Eig. 
1^ as a function of and AM, for 3000 fb“^ of inte¬ 
grated luminosity at LHG 14. The blue regions show the 
3(7 reach, while the red regions show the 5(7 reach. The 
highest reach is obtained for A M ^ 80 GeV. Eor larger 
A M, the significance drops due to reduced branching to 
leptons, as the dominant decay mode reverts back to the 
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FIG. 6: The left panel shows the pt distribution (normalized to unity) of the leading and next-to-leading leptons for 
signal (red unshaded histogram) and background (grey shaded histogram) after preselection cuts. The right panel 
shows the distribution (normalized to unity) of for signal (red unshaded histogram) and background (grey 

shaded histogram) after lepton selections. All masses are in GeV. 




FIG. 7: The left panel shows the pt distribution (normalized to unity) of the leading 6-jet for signal (red unshaded 
histogram) and background (grey shaded histogram) after the selection. The right panel shows the 

distribution (normalized to unity) mT 2 (right panel) for signal (red unshaded histogram) and background (grey 
shaded histogram), after all other selections. All masses are in GeV. 


usual t hW'^. For smaller AM, the 6-jets coming 
from stop decay become too soft for the signal events to 
pass the preselection cuts. 


B. Comments on Systematics 

The mass reach for our study shown in Fig. does 
not take into account the effect of systematics. We now 
make some comments in this regard. 

• Glearly, one of the most crucial factors determining 
signal discrimination is the identification of 6-jets, 
as part of the preselection cuts. This requires effi¬ 
cient and robust 6-tagging at low pr- In our study, 
we have taken the conservative 6-jet pt threshold 
of 30 GeV. The challenge of 6-tagging at low pr 
is expected to be more difficult in the pileup con¬ 
ditions of 14 TeV, although preliminary detector 
upgrade studies have shown the ability to go down 
to Pt = 30 GeV [46] . In the 8 TeV run, both GMS 


and ATLAS have shown the ability to identify 6-jets 
down to Pt = 20 GeV US]; this would improve the 
mass reach for compressed stop searches substan¬ 
tially. 

• Systematic uncertainties are expected to lead to 
substantial degradation of the mass reach. The 
background estimation and signal extraction would 
largely depend on the upgraded detector configura¬ 
tion and trigger conditions at 3000 fb“^, and the 
final systematics in the high pile-up environment 
would depend on the ability to reject pile-up. 

• The distributions in Fig. and Fig. show that 
significant improvement can be achieved by pursu¬ 
ing a shape based analysis. Although we do not 
perform such an analysis here, we give some esti¬ 
mates. Our S/B for the benchmark point in Table 
(njis similar to the S/B in sl compressed stop study 
performed by some of the current authors [39] , Ta¬ 
ble II. In that study, a shape based analysis of the 
pY distribution was performed using a binned likeli- 
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hood following the test statistic based on the profile 
likelihood ratio, with systematic uncertainties be¬ 
ing incorporated via nuisance parameters following 
the frequentist approach. Based on that study, we 
can estimate that the combination of a shape based 
study and incorporation of 3% systematics would 
lead to a reduction of ^ 0(100) GeV in the mass 
reach in Fig. which was conducted without the 
shape analysis and without systematics. Thus, we 
can expect a 5cr reach of around ^ 240 GeV, 
for AM ^ 80 GeV. We leave a more detailed study 
incorporating shape analysis for the future. 



FIG. 8: The significance where S and B are 

the signal and background rates, respectively, as a 
function of and AM, for 3000 fb“^ of integrated 
luminosity at LHG14. The blue contour shows the 3cr 
reach, while the red contour shows the 5cr reach. All 
masses are in GeV. 


VII. CONCLUSIONS 

The anomalous magnetic moment of the muon is one 
of the most precisely measured quantities in particle 
physics. The upcoming New g — 2 Experiment at Fer- 
milab will measure it with even greater accuracy, further 
establishing the discrepancy between the SM prediction 
and experimental value. Simple but powerful order-of- 
magnitude calculations of the effects of new physics show 
that supersymmetry is a prime contender in the resolu¬ 
tion of this discrepancy; the MSSM indeed provides am¬ 
ple parameter space in this regard. 


We have studied the regions of supersymmetric param¬ 
eter space preferred by the muon g — 2 anomaly in the 
context of constraints arising from the first run of the 
LHC. We have found that the main constraints arise from 
slepton searches in the dilepton channel, and chargino 
searches in the trilepton channel. A combination of these 
searches rules out a substantial part of the parameter 
space for tan/3 ^ 0(10). Most of the parameter space 
at larger tan^d ^ 0(40) remains unexplored, and the 
constraints here will depend on the future reach of slep¬ 
ton and chargino searches in the high-luminosity LHC. 

Since a primary motivation for supersymmetric 
searches at the LHC is the naturalness paradigm, we 
have evaluated the allowed parameter space for the muon 
g — 2 anomaly in the light of stop-assisted chargino 
and neutralino production. We have found that current 
constraints coming from stop searches in the i ty? 
and i hxf are already probing the entire parameter 
space at tan/3 ^ 0(10), and substantial parts of it at 
tan/3 ^ 0(40). 

The future performance of the LHC in these stop 
searches will be crucial in illuminating the interplay be¬ 
tween naturalness and a supersymmetric resolution to 
the muon g — 2 anomaly. One possible blind spot are the 
compressed stop scenarios. We have studied one such 
scenario, which can be dominant in the presence of light 
charginos and sleptons in the parameter space compati¬ 
ble with the muon 2 anomaly: i hiu{pi) hiuxi^ 
where AM = — m^o ^ 80 GeV. Using a simple cut 

and count approach, we have found that the 5cr mass 
reach in such a scenario is ^ 350 GeV. Systematics 
will degrade this reach substantially; to obtain sensitivity, 
we have pointed out that shape analysis of kinematic dis¬ 
tributions will be essential to discriminate between signal 
and background. This will further probe the regions of 
paameter space where naturalness and a supersymmetric 
resolution to the muon g — 2 anomaly are compatible. 
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